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Counting dark matter particles in LHC events 
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We suggest trying to count the number of invisible particles produced in missing energy events at 
the LHC, arguing that multiple production of such particles provides evidence that they constitute 
stable Dark Matter and that counting them could yield further insights into the nature of Dark 
Matter. We propose a method to count invisible particles, based on fitting the shapes of certain 
transverse- or invariant-mass distributions, discuss various effects that may affect the measurement, 
and simulate the use of the method to count neutrinos in Standard Model processes and Dark Matter 
candidates in new physics processes. 
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I. INTRODUCTION 

What a fillip it would be if the Dark Matter that 
abounds in the heavens could be manufactured here on 
Earth, at the LHC. A Dark Matter (DM) particle, being 
colour singlet and electrically neutral, would be invisible 



in the LHC detectors, but would manifest itself in the 
form of "missing energy" , recoiling against visible mat- 
ter. 

While production of DM would imply a missing energy 
signal, the converse is not obviously true. For example, 
an invisible particle produced at the LHC need only live 
long enough to escape the detectors, while the DM in 
the cosmos has been around for billions of years. The 
question would thus arise of how to establish the link 
between a new, invisible particle produced at the LHC 
and the DM in the cosmos. 

One way would be to measure the properties of the 
invisible particle (and any new companions) and thence 
to compute the relic density and compare with cosmo- 
logical observations. Alas, this is unlikely to be feasible 
at the LHC: even if we neglect the difficulties associated 
with making precision measurements at hadron colliders, 
there remains a basic obstacle, in that while the relic 
density is determined (presumably) by the weak interac- 
tions of new particles, production of those same particles 
at hadron colliders proceeds dominantly via the strong 
interaction. 

Here we should like to make a more modest, but per- 
haps more feasible, proposal for strengthening the hy- 
pothesis that a new, invisible particle produced at the 
LHC really is DM. Our proposal is simply to count the 
number of invisible particles in missing energy events. To 
begin with, our system of counting will be loosely based 
on the "one-two-many" system of the Amazonian Piraha 
tribe but simplified to "one- many". That is, we pro- 
pose to try to establish that invisible particles are being 
multiply produced in events. As we argue in the next 
Section, this constitutes evidence for the existence of a 
symmetry that stabilizes DM. 

If we can establish that new invisible particles are be- 
ing multiply produced, then it makes sense to follow the 
lead of the Piraha and to try to extend our system of 
counting beyond "one-many" . For example, we might 
try to count modulo some integer. We shall sec, for ex- 
ample, that even counting only modulo two may lead to 
further insights about the nature of DM particles (and 
antiparticlcs) and the symmetry that stabilizes them. 

Having argued, in S^jlll and lllll that counting invisible 
particles (or at least establishing that they are multiply 
produced) is a worthwhile thing to do at the LHC, we 
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sketch, beginning in ^IVl a method by which one might 
hope to do it. In a nutsheh, our strategy is to identify 
observables that depend strongly on the number of invis- 
ible particles present, but which are relatively insensitive 
to all the other variables over which, a priori, we have 
little control. We begin by considering, for simplicity, the 
limiting case in which the masses of the invisible parti- 
cles are negligible on collider scales. In ijIV Al we dis- 
cuss the case in which new coloured particles are singly 
produced at the LHC and subsequently decay to visible 
and invisible particles; in ijlVBi we discuss the case of 
pair production. In we discuss various complications 
that arise in the massive case. In WI[ we discuss a va- 
riety of other effects, including matrix elements, finite 
widths, backgrounds and upstream transverse momen- 
tum. In Win wc apply our method (neglecting show- 
ering and detector effects) to three Standard Model test 
cases, namely decays of a ly-boson to a charged lepton 
and a single neutrino, decays of the Higgs boson to a pair 
of charged leptons and a pair of neutrinos via interme- 
diate W-bosons, and decays of a pair of top quarks in 
the di-leptonic channel, and to a supersymmetric decay 
chain. Section IVIIII provides a summary of the main re- 
sults, while a number of the more cumbersome formulae 
have been relegated to an Appendix. 



II. MULTIPLE PRODUCTION AND 
STABILIZING SYMMETRIES 

Why is multiple production of invisible particles ev- 
idence for DM? The obvious way to make DM suffi- 
ciently long-lived is to stabilize it by means of a sym- 
metry: If DM is the lightest particle transforming non- 
trivially under the action of the symmetry (henceforth 
"non-singlet"), then it cannot decay. ^ Suppose (as seems 
likely) that the dark matter particle is heavier than the 
constituent quarks of the proton, such that the quarks 
are singlets. Then protons and pairs of protons are sin- 
glet states, as arc the products of LHC pp collisions. If wc 
arc lucky enough to produce a non-singlct DM particle, 
then such a final state must contain another non-singlct 
particle. In the simplest case, this would be a second 
DM particle (or its anti-particle), but it could also be a 
different particle. If it is different, and if it is both visible 
and stable on detector length scales, then we will not see 
multiple production of invisible particles, but we will see 
a spectacular charged track. In all other cases, we must 
end up with multiple invisible particles in the detector, 
in order to form a singlet final state. 

Before going on, let us argue that the observation of 
multiply-produced invisible particles would not be a triv- 



ial result, in that the alternative hypothesis of a new, 
singly-produced, invisible particle can plausibly be enter- 
tained.^ Indeed, one might worry that if one could draw 
a Feynman diagram leading to significant production of 
such a state in association with SM states at the LHC, 
then the crossed diagram with only the invisible particle 
in the initial state would imply a lifetime less than the 
time required to traverse an LHC detector. However, this 
is not necessarily the case. 

Consider for example a resonance with mass M, 
strongly produced at the LHC, which decays into an 
electrically-neutral and colour- neutral particle of mass m. 
The crossed Feynman diagram, with virtual exchange of 
the resonance, can contribute to the decay of the "in- 
visible" particle through a process which, at worst, in- 
volves a three-body final state. Then the lifetime is 
T ^ 2567r^M^/m^. Being produced with typical mo- 
mentum of order M, the "invisible" particle will travel 
an average distance L ^ rM/m before decaying. A pre- 
dominantly invisible decay (say with i > 30 metres) is 
obtained for m < 2 GeV(Af/TeV)t. This allows for a 
reasonable range of masses. Note, in particular, that 
the particle can be sufficiently heavy to avoid constraints 
from stellar cooling. Moreover, the bound on m gets 
weaker if the production process occurs in association 
with heavy SM particles (such as W, Z, t, or 6), since 
the decay of the invisible particle must involve multi- 
body final states. 

III. COUNTING DARK MATTER 

If we can establish that invisible particles are being 
multiply produced, we may also gain further insights by 
counting them, perhaps modulo some integer. As an ex- 
ample, if we could count the invisible particles modulo 
two, then observation of an odd number of invisible parti- 
cles would rule out the simplest stabilizing symmetry, viz. 
Z2. Indeed, an odd number of charged particles cannot 
form a singlet of Z2. The simplest remaining candidate 
for a stabilizing symmetry would then be Z3. An al- 
ternative approach to distinguishing Z2 from other sym- 
metries was discussed in 0, Q.^ Moreover, if we make 
the plausible assumption that only the DM particle is 
stable on detector length-scales, then observation of an 
odd number of invisible particles also tells us that either 
the DM particle cannot be its own antiparticle, namely a 
real boson or a Majorana fermion, or that the stabilizing 
symmetry cannot be Abelian. 

Indeed, any Abelian group is a product of U{l)s and 
finite cyclic groups and DM must be charged under at 



There are obvious analogues among the known particles: the 
proton is stabilized, perhaps accidentally, by baryon number, the 
lightest neutrino by fermion number, and the electron by electric 
charge. 



^ Neutrinos give a known example of invisible particles that can be 
singly produced, but such events are easily picked out, assuming 
conservation of lepton number, by the presence of net visible 
lepton number in the final state. 

^ For models with stabilizing symmetries other than Z2, see 
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least one of these. If, on the one hand, DM were its own 
antiparticle and were charged under a C/(l), one could 
not write the necessary mass term for DM in the La- 
grangian. If, on the other hand, DM were charged only 
under some Zjv, one could write the mass term but could 
not simultaneously form a singlet final state from an odd 
number of DM particles. 

We note, moreover, that counting modulo two is unaf- 
fected by the presence of neutrinos in the final state, since 
these too can be counted modulo two, via the visible lep- 
ton number in the final state. If we try to be even more 
ambitious and count invisible particles modulo N > 2, 
then we should need to solve the problem of counting 
neutrinos on an event- by- event basis. One could imagine 
doing so, by using our knowledge of neutrino dynamics, 
but we shall not go into the details here. 



IV. LIGHT INVISIBLE PARTICLES 
A. Single production 

Having argued that counting invisible particles (or at 
least establishing that they arc multiply produced) is a 
worthwhile thing to do at the LHC, let us now try to con- 
vince the reader that it can be done, at least in principle. 
To do so, we recall that invisible particles can be discov- 
ered by their recoil against visible particles, manifesting 
themselves in the form of missing energy in an event. 
Now, the total collision energy is shared out between all 
the visible and invisible particles in an event in a random 
fashion; as a result, the shape of the distribution of the 
missing energy in a sample of many events will depend 
on the number of invisible particles present. 

In fact, the distribution of almost any observable will 
depend on the number of invisible particles present. This 
is, of course, good news. The bad news is that almost any 
observable will also depend on many other things, making 
it hard for one to be sure that one is really counting the 
number of invisible particles and not just measuring some 
poorly-defined combination of many other things. This 
leads us onto the important issue of how we propose to 
count invisible particles in practice. 

To do so, one would like to find an observable which 
depends strongly (and in a known way) on the number 
of invisible particles, but which is less sensitive to all the 
other variables over which, a priori, we have little con- 
trol. Examples of such pernicious variables include: the 
sizes of Standard Model backgrounds, other unknown pa- 
rameters in the new physics Lagrangian (such as particle 
masses and couplings), the details of the event topology, 
the unknown boosts of produced new particles with re- 
spect to the laboratory frame, uncertainties in the parton 
distribution functions (pdfs) &c. There are other impor- 
tant effects to which we have less hope of being insen- 
sitive. These include: the unknown particle spins and 
widths, the details of hard process matrix elements, the 
presence of radiation in the initial state, detector reso- 



lution, and so on. Nevertheless, we shall argue that the 
sensitivity to these effects can be mitigated, for three rea- 
sons. One is that we can, of course, resort to numerical 
simulations of physics and detectors to try to model such 
effects. The second is that, as we shall see, the effect on 
certain distributions of changing the number of invisible 
particles is a dramatic one, and is unlikely to be masked 
or faked by other effects. The third reason is that, since 
we are trying to measure a discrete integer, we can toler- 
ate a large error and still have confidence in our result. 

Before proceeding to a more detailed discussion of 
these issues, let us begin by considering a simple exam- 
ple that illustrates the dramatic effect on distributions 
that may result from changing the number of invisible 
particles. We consider decays of a single parent parti- 
cle of mass M into some number of visible and invisible 
daughter particles. For the time being, we assume that 
the masses of the invisible particles are negligible at the 
scales of interest. We should like to define an observable 
which has a strong dependence on the number of invisible 
particles, which can be understood in a simple way. We 
now claim that (a version of) the transverse mass [To| 
is such an observable. We note that, given only a sin- 
gle visible particle, our observable could be any function 
of the momentum of the visible particle and the missing 
transverse momentum; with more visible particles, many 
more observables become possible, including, for exam- 
ple, the invariant mass of some of the visible particles. As 
described above, the distributions of any of these observ- 
ables will depend on the number of invisible particles. 
But the transverse mass (and invariant masses, which 
we consider in j|IV A2p are special among these because 
not only do their distributions have a maximal endpoint, 
which is easy enough to pick out in data, but also the con- 
dition for an event to be at the endpoint can be expressed 
in a simple way, in terms of the invisible particles. For 
example, in the case of the transverse mass, we shall see 
that necessary and sufficient conditions to be at the max- 
imum are that all invisibles be parallel and transversely 
directed in the rest frame of the decay.'' As we shall see, 
the transverse mass has the additional advantage that its 
behaviour near its maximum is independent of whether 
or not the decay involves intermediate particles that are 
on-shell. 

To see all this, let us group the daughter particles of 
the decay into a visible system of transverse momentum 
Pt and invariant mass my and an invisible system (of 
n particles) of transverse momentum and invariant 
mass TO/. The transverse mass is naively obtained by 
projecting the invariant mass of all daughters onto the 



* To avoid later confusion, we stress tliat the conditions change if 
the invisible particles have non-negligible mass. 
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observable transverse directions, namely 



'^T, naive 



my 



+ mj + 2 y {fl + m]){pl + ml) - ■ pr) ■ (1) 



However, the invariant mass of the invisible system, m/, 
is not observable and thus nor is TOT,naive- However, we 
can simply replace the true value of mj by some fixed 
value in the definition. We choose to replace it by its 
minimum value (which, since ?TiT,naive is a monotonically 
increasing function of m/, preserves the property rriT < 
M). This minimum value is given by the sum of the 
masses of the invisible particles, which we have taken, 
for now, to be negligible. We thus arrive at our final 
definition of the transverse mass, namely^ 



■PT 



(2) 



We now examine the conditions that events must satisfy, 
in order to be at the maximum endpoint of the tot dis- 
tribution. Since the invariant mass of the parent may be 
expressed as 



2{EvEi -f -pT- qvqi) 



(3) 



where E denotes the energy and q the longitudinal mo- 
mentum component; we have that 



ij - 2ev {\ffT ^ 

-2{EvEi - eyei - qvqi), (4) 



where we defined the transverse energies via = p^ 
rn? . Using the relation 



{EyEif - (eye/ -I- qvqif = [eyqi - qv^if 



(5) 



it is clear that EyEj > eyei + qvqi- But then the last 
three terms on the right hand side of ^ are each negative 
semi-definite, such that 



with equality iff. 



eyqi — qv^i = and mj = 0. 



(6) 



(7) 



The former condition is invariant under longitudinal 
Lorentz boosts and implies either that qy = qi = 
or that the configuration is a longitudinal boost thereof. 
The latter condition is invariant under all Lorentz boosts 
and implies that all the invisible particles' momenta are 
parallel (not anti-parallel). Thus, the condition to be at 



^ The transverse mass thus defined is, in fact, the natural vari- 
able to use in this case, in the sense that it captures all of the 
information that is available from kinematics alone 




FIG. 1; Normalized, phase space distribution of for n £ 
{1, . . . , 4} for the process M — >■ P + nX, i.e. the decay of 
a parent particle of mass M into a visible particle P and n 
invisible particles X. 



the maximum of the mx distribution is that all invisibles 
be parallel, and that there exist a longitudinal boost to 
a frame in which they be purely transverse. (The gener- 
alization to the case of massive invisibles will be carried 
out in gVl) 

We thus see that the conditions to be at the maxi- 
mum have a strong dependence on n (since one must 
line up all n invisibles), but are independent of (i) the 
number of visibles and their masses and (ii) longitudinal 
boosts of the parent with respect to the lab frame. We 
can guess the form of the tot distribution near its end- 
point, for configurations in which the parent is produced 
without any transverse momentum and decays according 
to phase space considerations alone, in the following way. 
The dependence coming from the condition that the total 
momentum of the invisibles be purely transverse is inde- 
pendent of n and may be evaluated at n = 1; it comes 
from the Jacobian that arises in changing variables from 
the decay direction (which is uniformly distributed) to 



-2) 

linear in mj and hence we get a factor of (^1 — 
In all, we find 



m^ and gives a factor of (^I — jfi'j ■ The n depen- 
dence comes from the distribution of the invariant mass 
of the invisible particles, which, as can be seen in (|A2|) 
goes as TOj^" Near the endpoint (at small TOj), m'^ is 



n+l 



dm^ 



M2 



(8) 



Wc reproduce the explicit calculation of the tot distribu- 
tion for the cases of one and two massless visible particles 
in the Appendix; the resulting expressions are given in 
(|A6IIA8|) and (|A13IIAI5|) and distributions for the first 
few values of n are illustrated in Figs. [T] and [21 We see 
that the dependence on the number of particles is, in 
some sense, the strongest possible: for n = 1, the distri- 
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FIG. 2: Normalized, phase-space distribution of for n G 
{1, . . . , 4} for the process M — ^ 2P + nX, i.e. the decay of a 
parent particle of mass M into two visible particles P and n 
invisible particles X. 

button diverges at the endpoint (in the absence of finite 
experimental resolution) , whilst it vanishes for n > 1 . So 
we can hope that it is possible to discriminate single and 
multiple production of invisible particles, by looking at 
the endpoint behaviour. (As regards the secondary ques- 
tion of counting the number of invisibles for n > 1, the 
fall-off with n is so rapid that a lack of statistics near the 
endpoint would be a concern for n greater than a few.) 

To sec whether such a hope could become a reality at 
the LHC, one must consider the various aforementioned 
complications. Happily, the result is robust with respect 
to a number of effects. We have already argued that, 
since the transverse mass depends on only transverse 
or Lorentz-invariant quantities, the phase-space distri- 
bution will be insensitive to pdfs (which give rise only to 
longitudinal boosts of the singly-produced parent) and 
their uncertainties. We have also argued that, at least in 
the limit of massless invisibles, the endpoint behaviour is 
independent of the number of visible particles and their 
masses.^ We now show that, in the same limit, the power- 
law fall-off in niT is also unchanged if we let intermedi- 
ate states be on-shell (leading to, for example, cascade 
decays) . 



1. Topology dependence 

It is easy enough to convince oneself that the endpoint 
behaviour of the tot distribution will be unaffected by 
the presence of on-shell intermediate states, with one 
caveat. Indeed, we have already seen that necessary and 
sufhcient conditions to be at the endpoint arc that the 



Even if there were a dependence, this would be easy to com- 
pensate for, since the number and masses of visible particles are 
observable. 




FIG. 3: Normalized, phase-space distribution of for the 
cascade decay process Af — > P -I- V — >■ 2P -\- X for various 
values of the intermediate mass ratio z = my /M^ . 

massless, invisible particles be parallel (to/ = 0) and 
transverse, qi = 0, regardless of the topology. The prob- 
lem then reduces to studying how easy it is to satisfy 
these two conditions, as one varies the topology. The 
condition that the particles be transverse cannot depend 
on the topology, since the phase space for any topology 
is always a Lorentz invariant. For the condition that the 
particles be parallel, an explicit calculation confirms that 
^ TOj^" as TOj — independently of the topol- 
ogy, with one exception. The exception is that there exist 
topologies for which one cannot hope to satisfy to/ — 0, 
namely those in which an on-shell intermediate decays 
exclusively into invisible particles. Then, the invisible 
particles from that decay cannot all be parallel, if energy 
and momentum are to be conserved.^ In all other cases, 
the endpoint behaviour of the tot distribution will be 
unaffected by the presence of on-shcU intermediates. 

In (jA25p . wc give the explicit form of the distribution 
in the simplest case of a decay involving two massless, 
visible particles and one massless, invisible particle, con- 
nected by a single intermediate state of mass my] one 
may easily show that the endpoint behaviour is given by 
([8]) for all my. Distributions for various intermediate 
mass ratios, z = ttiy/M^, are shown in Fig.[3l 



2. Invariant mass observables 

In cases where there are multiple visible particles 
present in the final state, one also has the option of us- 



^ On reflection, such a situation must result in a pathology, since 
even the most diligent experimentalist could not be privy to 
the knowledge that an invisible intermediate particle had de- 
cayed into other invisible particles. Clearly, what would happen 
in practice would be that the intermediate particle would be 
counted as a single (massive) invisible particle. 
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n=4 V^'*-..._^^^'-«.w=2 > 



FIG 



4: Normalized, phase-space distribution of iriy for n £ 
. , 4} for the process M 2P + nX. 



ing the invariant mass of some combination of the visible 
particles as an observable. Compared to the transverse 
mass, these have the advantage of being invariant under 
all boosts and we shall see that they too have endpoint 
behaviours that depend strongly on the number of invis- 
ible particles. Perhaps their most appealing property is 
that they do not require one to measure the missing en- 
ergy in events, with the well-known difficulties that such 
a measurement entails.^ 

Let us proceed to discuss how the endpoint behaviour 
of such observables depends on the number of invisible 
particles. Consider the case in which a single particle of 
mass M decays into / massless, visible particles and n 
massless, invisible particles. Now, it is simple to show 
that the invariant mass, my, of the I visible particles 
reaches its maximum only when all n invisible particles 
are produced with vanishing momentum in the rest-frame 
of the decaying particle.^ Thus, if we add one extra 
invisible particle, the extra condition to remain at the 
maximum of my is stronger than the extra condition to 
remain at the maximum of m^; to wit, in the first case 
the extra particle must have vanishing momentum (in 
the rest frame of the decay), whilst in the second case it 
need only be parallel with the existing invisibles. As a 
result the coefficient of n in the power-law fall-off near 
the endpoint is increased; in the simplest case / = 2, an 
explicit calculation shows that 



+2 



drriy 



1 - 



M2 



(9) 



We reproduce the full my distribution in (jA17p : distribu- 
tions for the first few values of n are illustrated in Fig. S) 



Invariant mass distributions have the apparent disad- 
vantage that, unlike the ttit distribution, their behaviour 
at their maximum docs change when intermediate par- 
ticles go on-shell and therefore use of them to count in- 
visible particles requires one to first know (or assume) 
the event topology. This is hardly surprising, since the 
location of the endpoint itself changes as intermediate 
particles go on shell. To give an explicit example, in 
the simplest cascade decay discussed above with two vis- 
ible particles and an on-shell intermediate of mass my, 
the endpoint is at my = — my and the distribution 
there is flat. But for a point-like three body decay, the 
maximum is at rriy = and the distribution vanishes 
linearly there. 



B. Pair production 

We now switch our attention to decays of pair- 
produced particles. Pair production is relevant in the- 
ories, such as minimal supersymmetric models, in which 
the parent states are charged under the stabilizing sym- 
metry and must themselves be multiply produced. It 
turns out that the mx distribution still has the same form 
of power-law fall-off for the case of pair- produced parents, 
assuming that the two parents are produced at rest, rel- 
ative to one another, and without an overall transverse 
boost. These are, perhaps, not such good assumptions at 
the LHC where, even in the absence of initial state radia- 
tion, two parents may be produced boosted back-to-back 
in their rest frame. One way to avoid this problem 
would be to use, instead of my, the invariant mass of the 
visible decay products coming from one of the decaying 
parents. This observable is, of course, invariant under 
any boost of that parent, but the endpoint behaviour 
of its distribution will only yield information about the 
number of invisible particles produced in that parent's 
decay. Alternatively, one may use the distribution of the 
observable tot2 [12, [3 ■ Again, this is the natural ob- 
servable for pair production, in that it encodes all of the 
information that is available from kinematic considera- 
tions [ll|, 13 More importantly as regards our 
discussion, it enjoys a large invariance with respect to 
boosts of the parent particles, reducing the sensitivity to 
the details of the production process. This invariance 
is exhibited most simply by noting that if we only con- 
sider events in which the invisible particles recoil against 
the visible particles in the decays, such that there is no 
initial state radiation or other source of transverse mo- 
mentum "upstream" in the event, then the algorithmic 



* Measurement of the missing energy might be required for trig- 
gering, or to establish the presence of a signal, however. 

® Provided the decaying particle is boosted with respect to the lab 
frame, there will still be missing energy on which to trigger, if 
required. 



Parents will necessarily be boosted if the production process in- 
volves non-zero angular momentum. 

Here, we focus on pair production of parents of identical mass; 
generalizations of mT2 might be more appropriate in other 
cases. 
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definition of 771^2 reduces to an algebraic expression |16l - 
[isj . Tlie expression depends on wfietlier a confi gura tion 
is "balanced" or "unbalanced" as defined in Il3| and 
is given by-'^^ 

2 j At + -y/^y — mfm^, balanced, 
I max(TO^, TO2), unbalanced, 

where At = \/v\t + "^w/ vir + ™2 +Pit •P2T andpi,2T 
and mi 2 denote the transverse momenta and invariant 
masses of the visible systems in each decay. 

This expression is manifestly invariant under inde- 
pendent longitudinal boosts of each parent, along with 
boosts that correspond to first boosting both decays to 
be purely transverse and subsequently applying equal, 
but opposite transverse boosts to each parent pTf."'^'^ 

Just like mx, the distribution of mr2 has an upper end- 
point encoding information about the masses of the par- 
ticles involved. However, the endpoint behaviour is com- 
plicated, if there are multiple visible particles in one or 
other decay, by the fact that either the unbalanced or bal- 
anced configuration may dominate, since both can reach 
the same maximum, M . For balanced configurations, the 
maximum is reached when the invisible particles are par- 
allel and transverse, leading to the same power-law be- 
haviour as in ([5]); for unbalanced configurations, we need 
to maximize one or other visible invariant mass, leading 
to the fall-off given in ^ , but with my — )■ •mT2 and with 
n replaced by the number of invisible particles in one or 
other decay. Now, since the endpoints are the same, the 
distribution there will be dominated by the component 
with the smaller power. For pair decays that are sym- 
metric in the sense of having the same number of invisible 
particles on each side, the balanced solution dominates; 
for asymmetric decays, the unbalanced solution always 
dominates (provided that the decay with fewer invisible 
particles has multiple visible particles). We note that 
in the asymmetric case, the power is always an integer, 
whilst it is strictly a half-integer in the symmetric case. 
Thus, no ambiguity can arise between the two. Just as 
for mT these results do not depend on the precise number 
of visible particles. 

In the simplest case of a single, massless, visible parti- 
cle in each decay, mTi reduces to 

•m%2 = \V\T\\V2T\^V\T -ViT, (10) 

which allows one to obtain an analytic form for the mT2 
distribution for two parents produced at rest in the lab or 
boosted as described above, given in (|A28p for the case 
of two parents of equal mass, each decaying to fc or /, 



As we did for (c/. ll2]|), wc have replaced the true, unknown 
value of the invariant mass of the invisible system in eaeli decay 
with zero in the definition. 

We stress that this invariance group, large though it is, docs not 
contain all back- to-back boosts of the parents. 




FIG. 5: Normalized, phase-space distribution of 771^2 for 
(fc,Z) e {(1, 1), (1,2), (2,2), (1,3)} for the process 2M 
2P+{k + l)X. 



massless, invisible daughter particles. The distributions 
for the first few values oi n = k + I are shown in Fig. [5l 



C. Counting the number of parents 

At this point, the reader might worry how a consci- 
entious experimentalist, blessed with a sample of sig- 
nal events, is to decide whether they came from de- 
cays of singly- or pair-produced parent particles and thus 
whether to try to count the invisible daughter particles 
using the uit or tot2 distribution. We now suggest two 
possible ways to proceed. 

The first possibility is simply to make an educated 
guess. For example, if the final state contains only a sin- 
gle visible particle, then single production seems likely; 
alternatively, if the final state contains two copies of all 
particles (or particle-antiparticle pairs), then pair pro- 
duction might reasonably be assumed. 

Perhaps the ideal strategy would be to try to indepen- 
dently count the number of parents and the number of 
invisible daughters. We have argued above that one ex- 
pects that the distributions computed here will only give 
a good fit when one chooses the "right" observable for 
the process. If one chooses the wrong observable, then 
the absence of the desired invariance properties will lead 
to an observed distribution which is rather poorly fit- 
ted by the idealized distributions we present here. One 
thus might hope that one can identify the right observ- 
able (and in doing so, count the parents) by comparing 
multiple observables and seeing for which of those one 
obtains an acceptable fit. As an example, for single pro- 
duction followed by decay to two visible particles and one 
invisible particle, one would expect a good fit to the tut 
distribution, but a poorer fit to the mT2 distribution. 
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V. HEAVY INVISIBLE PARTICLES 

Thus far we have assumed that the invisible particles 
have negligible masses; we now redress this. 

Changing to a non-zero mass has several ramifications. 
The first issue arises in the definition of the observables 
rriT and mT2- For massless particles, we saw that the 
natural definition was to replace m/ in ([1]) by its min- 
imal value, namely zero. For massive particles, the ob- 
vious analogue would be to replace mj by its minimum 
value, namely the sum of the masses of the invisible par- 
ticles present. Unfortunately, we do not expect that this 
quantity would have been measured at the time that the 
analysis proposed here could be carried out.^^ Again, the 
obvious solution, since jtlt as defined in ([T} is a mono- 
tonically increasing function of m/, is to replace to/ by 
its minimal possible value, namely zero. Thus we retain 
the definition This will lead, however, to yet further 
subtleties, as we shall see in the next Subsection. 

A second issue is that the phase space available, in 
the limit that the invariant mass of the invisible particles 
approaches its minimum, differs depending on whether 
those invisible particles are massive or not. This should 
hardly come as a surprise, since the condition to be at 
the minimum differs in the two cases: in the massless 
case, the momenta of invisible particles must be paral- 
lel, whereas in the massive case, the momenta should, 
in addition, vanish, in the CM frame of the invisibles. 
There is, thus, less phase space available near threshold 
in the massive case; the explicit behaviours are given in 
(|A2[) and (jA3p . As a result, ([5]) is changed in the massive 
case, for a point-like decay of a single parent particle into 
a single, massless, visible particle and n invisible parti- 
cles, to 



of parents, each into a single, massless, visible particle, 
wc find 



1 - 



"''T2 



(12) 



where now the cndpoint is given by M| 



TP' 



and n = k + l}^ Again, the behaviour is 
unchanged, compared to the massless invisible case, for 

These changes in the power law in going from massless 
to massive invisible particles may seem confusing when 
one considers that (fTTj). for example, holds for all non- 
zero values of /i„, however small. Naively, this suggests 
that there is a discontinuous (and observable) difference 
in the toj^ distributions at /i„ = and as /i„ — > 0. What 
really happens is that the behaviour in ([TT|) is obtained 
only in the region of the distribution where the invisible 
particle mass cannot be neglected in comparison with 
the distance from the endpoint. Further away, a tran- 
sition to the behaviour in (jS]) occurs. The location of 
the transition thus depends on the mass itself and the 
distributions show a smooth behaviour as one takes the 
limit in which the invisible masses vanish. We learn that 
there is no paradox, but we also learn that the endpoint 
behaviour we can expect to see is contingent upon the 
unknown mass of the DM. Our recommended strategy 
for measuring n will be to fit n in the endpoint region of 
the distribution, treating the unknown invisible mass (or 
masses) as a nuisance parameter in the fit. This is un- 
likely to result in a precise determination of the invisible 
particle mass, but it should give us enough flexibility to 
reliably fit n. 



+ 1 



dm^ 



1 - 



1/2 



(11) 



and where 



where now the endpoint is at Mi _ — jj- 

= 5]TOj5f is the total mass of the invisible particles. 
We note that the behaviour is unchanged compared to 
the massless case for n = 1. 

Unlike the case of massless invisibles, we do not ob- 
tain the same formula for the distribution for pair- 
produced particles. There, to be at the maximum of mT2, 
the invisible particles produced in one decay must be at 
rest relative to one another, but need only be aligned 
(transversely) with the invisible particles produced in the 
other decay. Consequently, for point-like decays of a pair 



A. Multiple visible particles 

A third issue in the case of massive invisibles is that, 
with mx defined as in the condition to be at the max- 
imum (and indeed its location and endpoint behaviour) 
changes, if there are multiple visible particles. To see 
this, let us revisit the arguments of i jIV Al Equation 
is unchanged, but since to/ cannot vanish, being bounded 
below by the sum of the masses of the invisible particles, 
necessary conditions for tot to be at its maximum are 
that eyqi — qv^i = 0, as before, and for m/ to be at its 
minimum, to/. Then, ([¥]) reduces to 



TO^ 



2e\ 



(13) 



It is, by now, clear that there is enough information in kinematics 
(or equivalently my or mj'2) to measure the DM mass at the 
LHC |19| in principle; it is also clear that the measurement will 
be extremely challenging in practice. We therefore take the view 
that the DM mass will be unknown at the time that efforts to 
count the number of invisible particles begin. 



Wc stress that this result is derived for parents produced at rest, 
or boosted as described after I llOt . If one includes arbitrary back- 
to-back boosts, the endpoint shifts to — min(/i(; , ) UM , but 
the change in the endpoint behaviour is observed in simulations 
to be not greatly changed. 
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FIG. 6: Phase-space distribution of mr for the process 2M — s> 
2P + IX, for a massive invisible particle of mass mx for var- 
ious values of the mass ratio x = m\/M'^. 



The expression in parentheses never vanishes and so to be 
at the maximum of tot one must also be at the minimum 
of ev and ergo my. Thus, to be at the maximum of to^, 
one must also line up all the visible particles (if they are 
massless), or make their momenta vanish in their CM 
frame (if they are massive).^® 

To illustrate this, we give in (|A23[) the to^ distribution 
for decay of a single parent to two massless visibles and a 
single massive invisible. The reader may verify that the 
power-law fall-off near the endpoint has changed from 
minus one-half (as in the case of a massless invisible) 
to plus one- half (as in the case of a massive invisible), 
because we now have to line up two (massless, visible) 
particles. Nevertheless, as illustrated in Fig.[6l the mass- 
less and massive distributions map onto each other as the 
invisible mass is sent to zero. 

Finally, we remark that this change in the endpoint 
behaviour, while dramatic, can easily be accommodated 
in performing the measurement of n, since the number of 
visible particles is easy to measure. 

Further subtleties arise for tot2, once we allow for mul- 
tiple visible particles. Most importantly, either the end- 
point at ^ , corresponding to balanced con- 
figurations, or the one at max(M — Hk,if , corresponding 
to unbalanced configurations, may be the true maximum 
in asymmetric decays. In the latter case, the endpoint 



This is one of two effects that enables the (sum of) invisible 
particle mass(es) to be measured: if, as is the case here, the 
hypothesized mass in the definition of my is smaller than the real 
mass, the maximum of niT is at the minimum of my , while if it is 
greater, the maximum of my is at the maximum of my , leading 
to a kink in the endpoint of the mx distribution as a function 
of the hypothesized mass 1201 . The other effect needs only one 
visible particle and relies on the varying alignment of the visible 
particle(s) with respect to upstream transverse momentum to 
generate a kink |21| . 



behaviour of the distribution is the same as that 

of the visible invariant mass distribution for the relevant 
decay, whose form wc now discuss. 



B. Invariant mass observables 

It is also worthwhile to consider the effect of massive 
invisible particles on the distributions of invariant masses 
of visible particles. Again, let us consider a single par- 
ticle of mass M decaying into / visible particles and n 
invisible particles, all with arbitrary (but non-vanishing) 
masses. Just as in the case considered before, where all 
the invisible particles were massless, one may show that 
the condition to be at the maximum of the invariant mass 
of the I visible particles is that all invisible particles be 
produced with zero momentum in the rest frame of the 
decaying particle. An explicit computation shows that 
the endpoint behaviour is given, for point-like decays, by 

m 

where now the endpoint is at A/3 = M — /Lt„. The coeffi- 
cient of n in the power-law fall-off is easy to understand. 
Here, unlike the massless case, if one adds an extra in- 
visible particle, then the extra condition to be satisfied 
by that particle to be at the maximum of either hit or 
TOy is the same: the extra particle must be at rest with 
respect to the invisible particles already present. Note 
however that the power laws are different in the massless 
and massive cases (2n — 1 and 3^ — 1 respectively), even 
though the condition to be at the maximum (all invisible 
particles' momenta vanishing) is the same in both cases. 



C. Topology dependence 

All the results obtained thus far for massive invisible 
particles were derived assuming point-like decays. In the 
limiting case of massless invisible particles, we saw that 
the endpoint behaviour of the transverse mass variable 
was independent of the decay topology. We now show, 
by means of a counter-example, that this is no longer true 
in the case of massive invisible particles. To wit, consider 
a single particle decaying into two invisible particles and 
a single visible particle. The behaviour of the transverse 
mass near its maximum is determined in part (as wc saw 
in W A[) by the behaviour of the invariant mass of the in- 
visible system near its minimum. For a point-like decay, 
the differential phase space in the square of the invisible 
invariant mass vanishes as the square root near its min- 
imum, whereas it is flat if an on-shell intermediate lies 
between the two invisible particles. Happily, we would 
still fit the correct number of invisible particles, namely 
two. if we incorrectly assumed a point-like topology. We 
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would, however, obtain an incorrect value of the mass, 
namely zero. 

Thus, when invisible particle masses are non- negligible, 
both transverse and invariant mass distributions may 
exhibit a dependence on the invisible particle mass (or 
masses) and on the decay topology. As such, we rec- 
ommend that, since neither of these is known a priori, 
they both be included as nuisance variables in the fit. 
Of course, one may hope that, in many cases, the de- 
pendence on these nuisance parameters will be small and 
will not result in an ambiguity in extracting the number 
of invisibles. We shall see instances of this in Will As 
an example of how this may arise, the fact that trans- 
verse mass distributions are independent of the topology 
in the massless limit suggests that the dependence on the 
topology will remain small as long as the masses are not 
too large. 



VI. OTHER EFFECTS 

A. Matrix element and spin dependence 

We now turn to other effects. The first of these is 
the modulation of phase space distributions by matrix- 
elements, due to spin or other effects. By fitting the phase 
space distribution, we are implicitly assuming that these 
effects may be neglected. In this respect, our approach is 
reminiscent of the on-shell effective theory (OSET) ap- 
proach [H [li- 
lt is easy to see, however, that there are cases in which 
these effects cannot be neglected for our purposes, in 
that they can. in principle, fake the effect of changing 
the number of invisible particles. Imagine, for example, 
that angular momentum considerations force two invisi- 
ble particles to be always produced with equal momenta 
(or parallel if massless). Then, it is clear that the re- 
sulting distributions will be indistinguishable from the 
distributions that would have resulted from a single in- 
visible particle (with mass given by the sum of the two 
individual masses). However, such effects appear only 
when decays occur at threshold and even then will be 
washed out by finite width effects. As an explicit exam- 
ple, in the decay h —J- 2W — >■ 212^ of a Higgs boson of 
mass 2mw, the two neutrinos are parallel in the limit of a 
narrow-width Higgs, but the rriT distribution still shows 
the n = 2 endpoint behaviour ^ in practice [25j . 

A related worry is that the inclusion of spin effects 
will lead to two distributions, corresponding to different 
values of n, that are indistinguishable near their end- 
points, if not elsewhere. Again, it is easy to construct 
an example where this could happen. Consider, follow- 
ing [2^, a scalar particle. A, that decays to an invisible 
fermion, ip, and an off-shell, Dirac fermion, which in 
turn decays to an invisible antifermion ?/;, and a visi- 
ble scalar, B. Let the relevant Lagrangian couplings be 
+ B-^PriIj + h.c, where 2Pl^r = 1 ± 75. Then, 
conservation of angular momentum forces the matrix el- 



ement to vanish when the two invisible particles are par- 
allel, which, we recall, is a necessary condition to be at 
the maximum of the distribution. As a result, the 
observed fall-off in the vtlt distribution will result in a 
measured value of n > 2, if we try to fit the endpoint 
behaviour to ([H). 

As is well known (see e.g. [131), similar effects may 
also arise if we consider invariant mass distributions of 
neighbouring visible particles in cascade decays. Such 
behaviour arises because the invariant mass is a function 
of the angle between the two visible particles, as mea- 
sured in the rest frame of the intermediate particle; a 
polarization of the intermediate particle imparted by the 
first decay can then result in an angular dependence of 
the products of the second decay. Our expectation is that 
such effects will be larger for invariant mass observables 
than for or tot2- Indeed, we have already seen in 
W Al that. in the case of massive invisible particles, one 
must minimize both m/ and my in order to be at the 
maximum of m^. Thus, in order for spin effects to pol- 
lute the endpoint behaviour, they must appear in both 
the visible and invisible particle systems and they must 
not compensate each other. 

Finally, it is possible that significant matrix-element 
effects will arise in another way. We have already seen 
that the transverse mass is maximized when the decay 
particles are produced perpendicular to the beam direc- 
tion. One can imagine that such configurations will be 
preferred or disfavoured by the hard production process. 
An obvious example occurs in Drell-Yan production of 
PF-boson, whose polarization disfavours larger values of 
rriT [13 (though the power-law fall-off is unchanged). 
However, we expect that significant effects of this kind 
are unlikely to arise at the LHC, where several hard pro- 
cesses typically contribute to a given final state. 

There is little more that can be done about these kinds 
of pathology at the general level, other than to be aware 
of their existence and to note that instances of them are 
rather rare (necessary conditions to generate effects in 
cascade decays, for example, are given in [1^ and are 
seen to be rather stringent). Moreover, even if such ef- 
fects are present, they are typically small. As such, we 
are inclined to take the view that one should be delighted 
rather than disappointed if one is forced to contend with 
such a pathology in practice, since it will also give us 
the opportunity to make inferences about the matrix ele- 
ments for processes involving new particles. The best way 
to proceed in such a case would seem to be a painstak- 
ing comparison of explicit Lagrangian hypotheses (with, 
e. g. differing spins, couplings and numbers of invisible 
particles) with data. 

Finally, another special case is the one in which invis- 
ible particles do not have a single mass value, but rather 
form a nearly continuous spectrum of masses. Obvious 
examples are theories with large extra dimensions and 
the fundamental gravity scale close to a TeV, or unparti- 
cle theories possessing a conformal sector. Both of these 
predict signals in e.g. monojet plus missing energy chan- 



11 



ncls, which we might mistakenly interpret as single pro- 
duction of a heavy coloured particle which then decays 
to a single visible particle plus n invisible particles. We 
have checked that it should be relatively easy to distin- 
guish these cases. An obvious difference is that the gravi- 
ton/unparticle distributions do not feature an endpoint 
and have an extremely rapid fall-off in tot. They are 
therefore unlikely to be confused. 



B. Finite width effects 

A decaying particle inevitably has a finite width, given 
by a Breit-Wigner distribution. This has two important 
effects. Firstly, the distributions above (derived for fixed 
mass and zero width) must be convolved with a Breit- 
Wigner distribution. Since the latter distribution is non- 
vanishing for arbitrarily large masses, the singular end- 
points discussed above disappear, to be replaced by dis- 
tributions which fall sharply beyond the peak mass. To 
illustrate this, consider convolving the phase-space mx 
distribution for a two-body decay (|A6[) with the Breit- 
Wigner distribution ((77 
suit is 



MY + rH'I^y'; the 



cos 9 



oc 



where tan 29 = 



(1- 



17 



'r ^2 



r2 
TP 



(15) 



A/2 



We plot the distribution for various values of the frac- 
tional width r/M, in Fig. [71 together with the corre- 
sponding distribution for two invisible particles. The 
presence of finite-width effects of this type means that 
it makes no sense to try to count invisible particles by 
directly fitting empirical distributions to phase space. 
Rather, we should fit the phase space distribution (with 
the normalization unfixed and the invisible mass kept as a 
free, nuisance parameter), convolved with a Breit-Wigner 
distribution (whose fractional width is also, in general, a 
nuisance parameter), in the region of the endpoint. 



C. Backgrounds 

It was argued in [3(| that the transverse mass ttit (and 
its cousin TOT2) also has the desirable property that it 
provides a good discrimination between signal and back- 
ground events, at least if the new physics signal is suffi- 




FIG. 7: Phase space distribution of mr for the process M 
P + nX, with fractional width y = V /M, for n = 1, 2. 



ciently heavy.^® To wit, whilst the signal can have large 
values of to^, background events have small values of 
mx- This property is not only useful for discovering new 
physics (and indeed the mT2 distribution in the dijcts 
plus missing energy channel using the 2010 LHC data 
gave the strongest constraint on squark masses [Tsl. [3ll|). 
but it also tells us that the upper endpoint of the tot 
or TOT2 distributions, where we propose to carry out the 
measurement of n, will have a suppressed background 
contamination. 



D. Upstream transverse momentum 

The Tot and tot2 variables (and, ergo, their distribu- 
tions) are both invariant under longitudinal boosts of the 
rest frame of the decaying parent or parents. The tot2 
distributions in (|A28|) are, furthermore, invariant under 
a larger set of boosts as described in §IVBI Neither vari- 
able is invariant under overall, transverse boosts of the 
(pair of) parcnt(s), however, and we now proceed to con- 
sider the effect of these. An additional benefit of the 
original definition of tot [13] as used to measure the mass 
of the W-boson [H,!!!] is that there is no shift in tot un- 
der a transverse boost /3t , at linear order in /3t HI] . Any 
shift thus goes as and is small for production of heavy 
states near threshold. What is more, since the endpoint 
of the Tot distribution is invariant under Lorentz boosts, 
any shift in the TOt distribution must also vanish as some 
power of the distance from the endpoint. There is thus 



The precise result depends on what one means by "the phase 
space TTiT distribution". Here we have used , ^ = , which 

corresponds to the decay rate under the assumption that all 
kinematic and dynamic dimensionful parameters are set by the 
variable parent mass and arc distributed according to the Breit- 
Wigner formula. 



More precisely, the mass difference between the parents and 
the invisible daughter in the new physics signal should be large 
compared to mass scales associated with Standard Model back- 
grounds. 

For decays involving jets in the final state, radiation from the ini- 
tial state poses an additional combinatorial ambiguity. Methods 
to reduce this using the 171^2 variable were discussed in [s^l . 
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a double suppression near to the endpoint. 

Unfortunately, neither of these results holds generally 
when we extend to massive final state particles (vis- 
ible or invisible). Indeed, one may easily show that the 
shift in rriT (as defined in ([T}) under a transverse boost 
vanishes at linear order only if the correct invisible mass 
is used in the definition of rriT and only if m/ = my. 
Similarly, it is easy to show that the upper endpoint of 
the rriT distribution is Lorentz invariant only if the cor- 
rect invisible mass is used in the definition of niT-'^^ Since 
we do not know the invisible mass a priori, this cannot 
be done. 

In conclusion, the distributions considered here will be 
sensitive to the presence of upstream transverse momen- 
tum in the event, in whatever form. Again, we hope the 
resulting distortions of distributions will not be so severe 
as to lead to a mismeasurement of the number of invisi- 
ble particles, but caution should be exercised. Of course, 
if one did have confidence in one's understanding of the 
spectrum of upstream transverse momentum associated 
with a new physics signal, then one could compensate for 
its effects. 



VII. EXAMPLES 

We now apply our fitting procedure to various ex- 
amples including both masslcss and massive invisibles. 
These test cases feature many of the effects discussed 
above and, as such, serve to verify the viability of our 
algorithm for counting the number of invisibles in the fi- 
nal state in the presence of, for instance, pdfs, cascade 
decays, finite widths and spin effects. 



A. Light invisible particles 

For realistic examples with massless invisibles, we need 
look no further than the Standard Model. Specifically, we 
consider M^-boson production and Icptonic decay at the 
Tevatron; Higgs boson production and di-leptonic decay 
via a pair of intermediate W-bosons at the LHC, with 
rrih = 180 GeV; and finally pair production of top quarks 
undergoing di-leptonic decay at the LHC. 

We generate events in each case at tree level using 
MadGraph4/MadEvent [s^. We analyze these at parton 
level, using truth information to avoid combinatorial am- 
biguities, e.g. in pairing the 6-jets with leptons in pairs 
of top quark decays. We then fit the transverse mass 
(or mT2) distributions in the region of the endpoint us- 
ing the relevant power-law behaviour for massless invis- 
ibles (summarized in Table HIT)) . convolved with a Breit- 









i?-squared 




Process 


Observable 


n = 1 


n = 2 


n — 3 


W -y If 


in-T 


0.982 


0.723 


0.529 


h -> W+W^ l+l-yv 


vaT 


0.501 


0.977 


0.779 




mrp2 


0.630 


0.996 


0.952 


ti W+W^bb i+l^yvbb 


mT2 


0.458 


0.999 


0.971 



TABLE I: Goodness-of-fit values for various Standard Model 
processes with neutrinos, fitted with the relevant endpoint 
power law (see Table llllf) . convolved with a Breit- Wigner. 
Fractional widths of parents are fixed to their true values. 
Best fits are shown in Fig. [S] 

Wigner distribution, whose fractional width is fixed to 
the true value for simplicity. The resulting i?-squared 
values, quantifying the goodness of fit, for various n are 
shown in Table U and the best and next-best fits are su- 
perimposed on the "data" in Fig. [8] 

We see that our simple algorithm picks out the cor- 
rect value of n in each instance, although the i?-squared 
values for n = 2 and n = 3 are sometimes close. Further- 
more, although the n = 2 distribution is the best fit for 
the rriT distribution in ft, — !• W~^W~, the quality of the 
fit is not quite as good as in the other cases and the ac- 
tual fall-off near the endpoint is somewhat steeper than 
n = 2 would warrant. We believe this to be a hallmark 
of spin correlations: in the limit that rrih — > 2mw and 
the widths go to zero, conservation of angular momen- 
tum forces the two neutrinos to be emitted in the same 
direction, hence acting like a single invisible. We see 
a similar effect in tt production if we pretend that the 
charged leptons, rather than the neutrinos, are invisible. 
This too is as expected, since spin correlations in the ti 
system have a stronger effect on the distribution of the 
charged leptons [36j . 

Finally, we note that for ti production we obtain a good 
fit using the endpoint power-law for point-like decays, 
even though the intermediate VK-bosons are produced on- 
shcU. 



B. Heavy invisible particles 

Processes with heavy invisibles in the final state 
are abundant in theories going beyond the Standard 
Model. We choose to fit an example where both 
spin correlations and topology-dependence are known 
to play a role, namely the supersymmctric process of 
ul pair production, followed by decays ul — ulX2 ~^ 
'^Lij/eR WL^A^flXi [IS- We use SOFTSUSY ^ to 
compute the SUSY spectrum, generate the events using 
Madgraph/Madevent and BRIDGE [111 and analyse them 



Again, it is this phenomenon that is responsible for the kink in 
the endpoint of the tut distribution that allows the masses to be 
measured with just one visible particle |21| . 



In fact we have verified that even for = 2mw , the best fit 
still corresponds to the power law for n = 2. 
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(a) m|, distribution for W production and leptonic decay at 



the Tevatron. 
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at the LHC, with = 180 GeV. 
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(d) rrtyj distribution for tt production and leptonic decay at 
the LHC. 



FIG. 8: Standard Model examples with final state neutrinos, along with the best and next-best fits for distribution of transverse 
mass variables near endpoint. The fits are to the relevant endpoint power law (see Table HIT)) convolved with a Breit-Wigner, 
with fractional width set to the true value. Details of Monte Carlo event generation are given in the text. 



using truth information as before. For the fits we use 
RAMBO to generate phase space for the point-like 
decays M + M ^ {3 + 3)P +{k + l)X with varying invis- 
ible mass; these are again convolved with a Breit-Wigner 
of fixed, true fractional width. The resulting i?-squared 
values for different n = k + l and invisible mass are shown 
in Table im with the best and next-best fits superimposed 
on the "data" in Fig. [HI Note that in spite of spin corre- 
lations that significantly affect the invariant mass distri- 
butions for pairs of visibles in this decay, once again the 
best fit for 171x2 is for the correct value of n, and also for 
the invisible mass that is closest to the input value of 115 
GeV. Moreover, as in the ti case, point-like phase space 
gives a very good fit to the distribution, even though the 
actual process is a 3-step cascade decay. Unfortunately 
it is not possible to say definitively, using the 771^2 distri- 
bution alone, that this process is not instead n = 3 with 
light invisibles. This is an ambiguity that comes about 
because the endpoint behaviour of 771,^2 here is governed 
by the unbalanced solution which falls off linearly; it can, 
however, be resolved by separately examining the visible 
invariant mass distributions on each leg of the decay (see 



i?.-squarod (Tn^2) i?-squarcd (^n^i) i?.-squarcd ("^v2) 



mx 1 GcV 
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fc = 2 
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0.947 


0.992 


0.899 


0.992 


0.855 


0.987 


0.843 


100 


0.998 


0.889 


0.822 


0.993 


0.921 


0.998 


0.911 


450 


0.859 


0.991 


0.870 


0.987 


0.911 


0.994 


0.900 



TABLE II: Goodness-of-fit values for the supersymmetric cas- 
cade decay in [23|, with neutralino mass of 115 GeV, fitted 
to the endpoint behaviour of generated (3 + fc) + (3 + i)-body 
phase space, where k + l = n, with varying invisible mass and 
convolved with a Breit Wigner of fixed fractional width. Best 
fits are shown in Fig. [S] 

FiginUb)).-^^ Finally, one might worry that the results of 
the fits might change if one allowed the fractional width 
in the Breit-Wigner distribution to float. We would argue 
that plausible values for these, at least for small n, can 
be chosen by eye from the distributions, and, moreover. 



As expected, the best fits for rriy are not as good as those for 
transverse masses, since we are fitting point-like phase space in 
spite of the known, large dependence of these distributions on 
topology. 
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TABLE III: Endpoint power-law fall-offs for transverse- and 
invariant-mass distributions, for massless {fi = 0) and mas- 
sive (/i 7^ 0) invisible particles. Valid for point-like decays 
(arbitrary decay topologies for mr in the /i = case) and 
only for decays with a single visible particle for tut or 771x2 
in the 7^ case. The formula for asymmetric pair decays 
applies only if there are multiple visible particles in the decay 
with fewer invisibles; if not, the formula for symmetric pair 
decays applies. The mv observable may be employed for any 
production process, provided the visible decay products of a 
single parent can be isolated. 

that scanning over these is unlikely to make a significant 
difference to the best fit. (For example it is difficult to 
see how a peaked n = 1 distribution could be mistaken 
for n = 2, even allowing for arbitrary fractional width.) 



VIII. DISCUSSION 

The existence of DM in the cosmos is now established 
beyond reasonable doubt; the prospect that we may im- 
minently manufacture and be able to study DM in the 
laboratory is a tremendously exciting one, for particle- 
and astro-physicists alike. Here we have argued that a 
useful preliminary way to study an invisible DM candi- 
date at the LHC would be to count the invisible particles 
produced in an event and we have outlined a procedure 
by which such a measurement might be achieved. Do- 
ing so would not only increase one's confidence that one 
really had made DM (if multiple production could be es- 
tablished), but could also yield deep insights about the 
nature of DM and its stabilizing symmetry. 

To this end, we have shown that certain obscrvablcs 
have a strong, power-law dependence, in their endpoint 
behaviour, on the number of invisible particles; our re- 
sults are summarized in Table IIIIl This endpoint be- 
haviour follows from simple kinematic considerations 
(namely that the invisible particles should be produced 
either parallel to each other or at rest) and is, to some 
extent, universal. The universality, together with the 
strength of the dependence, means that counting invisible 
particles using the endpoint behaviour should be robust 
with respect to a number of conceivable complications 
and gives us high hopes that one could reliably count the 
number of invisible particles in practice, given favourable 
conditions. 

We have focused on transverse mass variables (and 
their cousins, such as tot2) and invariant mass obscrv- 
ablcs. Both feature maximal endpoints with a strong 



dependence on the number of invisible particles. The in- 
variant mass observables have the advantages that they 
do not involve the missing energy, such that the empiri- 
cal distributions may be cleaner, and simulations suggest 
that they have reduced sensitivity to invisible particle 
masses in cascade decays. In contrast, transverse mass 
observables have the advantage that they are less sen- 
sitive to the details of the decay topology (indeed, ttit 
is independent of it, in the limit that invisible particle 
masses are negligible) and that we expect them to be 
less sensitive to spin effects. 

There are important respects in which our results are 
not universal, which we have taken pains to point out. 
The most important of these is the effect of the mass 
of invisible particles. This effect may be significant and 
may even lead to ambiguities in counting the number of 
particles. For example. Table IIIIl shows that, for single 
production, one cannot distinguish (n = 3, = 0) from 
(n = 4, /i 7^ 0) using the power-law fall-off of my alone. 
There is also a secondary dependence on the widths of 
the parents involved and, in general, on the topologies via 
which they decay. To compensate for this, we propose 
that experiments allow these variables to float (within 
reasonable ranges) when fitting the endpoint behaviour. 
We have argued, nevertheless, that there will often be a 
high degree of insensitivity, in which case the poor con- 
straint on the value of the nuisance parameter that results 
may provide a useful cross-check. We have also identified 
pathological cases, in which matrix-element effects may 
lead to ambiguities or errors in inferring the number of 
invisible particles. 

In the massless invisible case, we have argued that it 
suffices to fit the endpoint power-law behaviour (con- 
volved with finite-width and, eventually, detector ef- 
fects). In the massive case, there is a transition between 
massive and massless power-law behaviours and so the 
approach needs to be modified. Concretely, our pro- 
posed strategy is as follows. Given some missing energy 
signal in a particular channel and associated transverse- 
or invariant-mass distributions, generate phase-space dis- 
tributions, corresponding to single or pair production of 
parents at rest in the lab frame, with the observed num- 
ber of visible particles. This should be done for varying 
numbers of invisible particles and their masses (and, if 
necessary, for various decay topologies). Convolve these 
distributions with Breit-Wigner distributions (of vary- 
ing width if necessary) and the appropriate detector re- 
sponse. Fit the resulting distributions to the signal in 
the endpoint region and extract the best fit value of the 
number of invisible particles. 

Our proposal might be further developed in various 
ways. Firstly, we have focused our attention on only 
a limited set of observables and only on their maxima. 
More generally, many observables feature singular points 
in their distributions that result from the projection of 
phase space into its observable subspace [40]; the be- 
haviour near these generalized singularities should also 
contain information about the number of invisible par- 
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at the LHC. Similar results arc obtained for the other leg. 

FIG. 9: Best and next-best fits for the SUSY cascade process in [23] with 6 visible particles and a pair of 115 GeV neutralinos 
in the final state. The fits are to the endpoint behaviour of generated (3 + fc) + (3 + i)-body phase space, where k + I = n, 
with varying invisible mass and convolved with a Breit Wigner of fixed fractional width. Details of the Monte Carlo event 
generation are given in the text. 



tides and merits further exploration. Secondly, we have 
not dealt with the issue of combinatoric ambiguities aris- 
ing in identical pair decays or from the presence of ini- 
tial state radiation. Variables generalizing 771^2 that deal 
with these (and which enjoy similar boundedness prop- 
erties) have been proposed [H, [s^l and it would be of 
interest to extend the study carried out here to them. 
Thirdly, we have shown that there can be a dependence 
on the topology in various cases and a fuller study of 
such effects would be desirable. Finally, even though we 
have yet to sec evidence for new, invisible particles pro- 
duced at the LHC; now would seem to be the ideal time 
for experiments to validate and refine our proposal, by 
counting the neutrinos which certainly have been abun- 
dantly produced in various SM processes, of which we 
have discussed various examples. 
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Appendix A: Useful Formulae 

The differential phase space for the decay of a particle 
M into n final-state particles Pi is 



{2tt 




Here P, E and pi, Ei are the momenta and energies of 
the initial and final state particles, respectively. In the 
case in which all particles Pi are massless, the integrated 
phase space is 



$„(Af ^ Pi...P„ 



l)!(n-2)!24"-5 7,2«-3■ 
(A2) 

In the case in which all particles Pi are massive, with 
masses m^, the integrated phase space is given, in general, 
by an elliptic integral. For the endpoint behaviour of 
distributions, we need only its behaviour near threshold, 
AI — > ^„ = YiiTTii, which is given by 



Pi 
1 



.P„ 



"n,mi 



2 2 TT 



^r(^) 



(A3) 



To compute the mx distribution for a single parent 
of mass M decaying at rest to a single, massless, visible 
particle (P) and n, massless, invisible particles (X), we 
note that the definition of mx reduces, in this case, to 



(A4) 



The differential phase space for the decay process M 
P + nX can be conveniently written by taking the con- 
volution of the two-body decay process with the system 
of n invisible particles (collectively denoted as /) , having 
invariant mass mj 







^ p- 



-nX) = 

d^2{M PI) d<i>n{I -> nX). (A5) 



' Using (jA2|) . we obtain the normalized phase-space dis- 
(Al) 

tribution in m^, in the center-of-mass frame (or a frame 
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boosted orthogonally to the transverse plane) 



such that 



where 
and 

An{z) 

such that 



1 



n(n- 1) /-i-^ . 



(A7) 



, for n > 1, 



A2(z) = log 



1 + VT 



(A8) 
(A9) 



i?2(w) = 2\/l - li; - (1 - w) + 4w log 



1- 



www 



1 — / IC' 

-(1 — 4:w)\l arctan 



1 — w 



jc(A16) 



The normalized phase-space distribution in the invari- 
ant mass of the two visible particles is 



$("+2) dm^ 



where 



(M^2P + nX) = C„[j^], (A17) 



A, (z) = 3 log ( ^ + _ 3x/r^, (AlO) 



Ci(a) = 2(1 - a) 



(A18) 



A4(z) =3(2 + 2) log ( ^ ) - 9VT- 



and 



(All) 

To compute the distribution with two visible particles, 
it is convenient to write the differential phase space as a 
convolution between the visible {V = 2P) and invisible 
(/ = nX) systems, having invariant masses my and to/, 
respectively. 



d$„+2(A/ -^2P + n 



dmydrrij 



<M (27r)2 



d$2(A/ ^ VI) d<p2iV 2P) d^n[I nX). (A12) 

The normalized phase-space distribution of the trans- 
verse mass ([2]) of the two visible particles, in the ccnter- 
of-mass frame (or boosted orthogonally to the transverse 
directions), is then 



M2 d$ 



n+2 



<^„^-2 dmi 



i'n+2 

where 

Bi{w) 
and 



{M -^2P + nX) = B„ ( , (A13) 



(1 - 2w) 
1 - arctan 

y/w{l - W) 
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(A14) 



(l-V^)(l-!//v^) 



Bn{w) =n^{n'^ -1) / dy 

Jo Jo 

a;"~2(l - y'^w-'^) 
2a/(1 - a;)2 - w - 2xy + y'^{l - w;-i) 



dx 

, for n > 1, 



(A15) 



Cnia) 



\t?-1) I dx x""^ a/(1 + a - a;)2 - 4a for n > 1, 

Jo 

(A19) 



such that 



C2(a) = 6(1 - a^) + 12aloga, 



(A20) 



C3(a) = 12(l-a)(l + 10a + a2) + 72a(l + a)loga, (A21) 



C4(a) = 20(l-a^)(l+28aW)+240a(l+3a+a^)loga, &c. 

(A22) 

For a point-like decay of a single particle to two, mass- 
less, visible particles and an invisible particle of mass 
mx , the TOy distribution is given by 



where 



17 



g {w, a;) + (1 — ui — 'Sx) ^^(1 — x)"^ — 2w {1 + x) + w 
D{w,x)^ ^ \ (i-^-^-)(i-^-^(3+2.-2..)) ^^^^^^ f^^±^\ forV^<l-^, 



and 



y/w{l-w} 



g X = 3xa/ (1 - x) - w + = ^ arccos . (A24) 

^ Vi«(l-u.) \ l^x-wj 



For a one-step cascade decay with two, massless, visible 
particles and a massless, invisible particle, the distribu- 
tion of rriT is given by 



/ ™2 ^2 



2 "'"''T 



(A25) 



where 



2(1 ~z)^VT^ 



dx- 



x{2^-x) 



{l + x^-y)^l - {^-xf 



(A26) 



For events in which two parents, both of mass M, are 
produced at rest and each decays into a visible particle 
[P] and to fc or / invisible particles, where k + l = n, the 
normalized phase-space distribution in is given by 



M2 di'^k+i^i+i) 



«>fc+i$/+i dm; 



T2 



where 



(2Af ^ 2P+{k+l)X) = F, 



^T2 



(A28) 



and where 




if y > 1 — z. 



(A27) 



= - [ da [ db 



and the integrals An{z) were defined above. 
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